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Abstract. We study the occurrence, visibility, and curvature strength of singularities 
in dust-containing Szekeres spacetimes (which possess no Killing vectors) with a 
positive cosmological constant. We find that such singularities can be locally naked, 
Tipler strong, and develop from a non-zero-measure set of regular initial data. When 
examined along timelike geodesies, the singularity's curvature strength is found to be 
independent of the initial data. 
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1. Introduction 

It has recently been shown that a central, locally naked, Tipler strong singularity 
develops in the gravitational collapse of inhomogeneous spherical dust with a positive 
cosmological constant |1|, |2|. Those results are a generalization to asymptotically 
de Sitter spacetimes of the well-known asymptotically flat spherical dust — Lemaitre- 
Tolman-Bondi (LTB) — collapse models. The exact solvability of such models enables a 
detailed analysis of gravitational collapse, and formation and structure of singularities, 
which shows that, from regular initial data, a Tipler strong curvature singularity 
always develops, which can be locally or globally naked, depending on the initial data 
0. Although almost a paradigm for gravitational collapse, LTB models rely on the 
simplifying assumptions of irrotacional dust and spherical symmetry. Within spherical 
symmetry, various generalizations of the matter model have been considered [§, ||, f| [7|,@], 
with similar results: a central curvature singularity forms, and its visibility (local vs. 
global) depends on the initial data. 

Comparatively, much less is known about singularity formation and structure in 
non-spherical collapse. Thome's seminal analysis of cylindrical collapse |J led him to 
formulate the hoop conjecture, which essentially states that horizons form if and only 
when the gravitational mass of the system is confined to a maximum "radius" in every 
direction [|10j. Subsequent numerical analyses of prolate and oblate collapse by Shapiro 
and Teukolsky [TTJ, and of gravitational radiation emission in aspherical collapse by 
Nakamura, Shibata, and Nakao pifl , were unable to refute Thome's conjecture. A 
perturbative analysis of LTB collapse by Hirada, Iguchi, and Nakao [O, [TR [R3], showed 



that linear non-spherical matter, metric, and matter coupled to metric perturbations 
remain bounded in the limit of approach to the Cauchy horizon. The very similar 
structure of the central singularity in LTB and Tolman-Bondi-de Sitter (TBdS) models 
0, [HJ, strongly suggests that analogous marginal stability properties also hold for the 



latter case. The quasispherical dust collapse models, given by the Szekeres metric [18 



were analyzed by Szekeres himself |19|| , Joshi and Krolak |2(| , and by Deshingkar, Joshi 



and Jhingan [pl[. As in the spherical LTB "central" singularity forms, which is 

Tipler strong, and can be locally or globally naked, depending on the initial data. 

Recent measurements of type la supernovae p2|, |23[| and peculiar motion of low 



redshift galaxies [f24|, appear to suggest the existence of a positive cosmological constant. 
The possibility of non-zero asymptotic (constant) curvature constitutes an obvious 
motivation for the study of gravitational collapse with a cosmological constant |2|, |25||. 
In this paper, we address one aspect of this problem by studying the formation, visibility, 
and curvature strength of singularities in quasispherical inhomogeneous spacetimes 
with a positive cosmological constant, thereby generalizing the results of [I]] to the 
quasispherical case, and those of (TP], |2T| to the asymptotically de Sitter case. 

The general metric for Szekeres spacetimes with a positive cosmological constant 
was obtained in closed form by Barrow and Stein-Schabes [p6fl , in the context of the 
"cosmic no hair theorem", but this solution also provides a valuable test-bed metric 
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for the analysis of singularity formation and structure in non-spherical gravitational 
collapse. One must remark, however, that such a departure from spherical symmetry, far 
from being arbitrary, is very well-defined — whilst the four- dimensional spacetime metric 
does not admit any Killing vectors, it does possess an invariant family of spherical two- 
surfaces, hence the name quasispherical. Accordingly, the results presented here provide 
a useful, but somewhat limited, insight into non-perturbative departures from spherical 
symmetry. Note also that, Szekeres spacetimes can be matched to a Schwarzschild (or 
Schwarzschild-de Sitter, if A > 0) spacetime, and thus cannot contain gravitational 
waves 



27 



As we shall see below, as in the Szekeres and TBdS cases, for quasispherical 
asymptotically de Sitter dust collapse, from regular initial data a central singularity 
develops, which is Tipler strong, and locally visible (for a large class of initial data). 
Our results suggest that neither lack of asymptotic flatness, nor "mild" deviations (in 
a sense to be precisely defined below) from spherical symmetry, play an important role 
in the formation and nature of singularities in gravitational collapse. 

The paper is organized as follows: Szekeres spacetimes with a cosmological constant 
are outlined in section 2. Section 3 discusses the existence and visibility of the 
singularity, along null and timelike directions. The curvature strength of the singularity 
is analyzed in section 4. Section 5 concludes with a discussion and summary. 

Geometrized units, in which G — c — 1, are used throughout. 

2. Szekeres spacetimes with a cosmological constant 

In this section we present the relevant equations for Szekeres metrics ]T8|, [HI with a 



positive cosmological constant [P6fl . The stress-energy tensor is 

Tab = ip + p)u a u h - (p - -^j g ab , u a = 5 a , (1) 

where u a , p and p are the four- velocity, pressure and fluid density, respectively. The 
metric can be written in normal Gaussian coordinates: 

ds 2 = -dt 2 + e 2a dr 2 + R 2 e 2u (dx 2 + dy 2 ), (2) 

where a = a(t, r,x,y), R = R(t, r), v = u{r, x, y), and 

e <*(t,r,x,y) = r ) + R ^ r y^ X; y y + fc( r )]-£ ; ( 3 ) 

e -u(r^y) = Ci{x 2 + y2) + 2(caX + C3y) + (4) 

k(r) = A(cxC A -c\- c\) - 1, (5) 

where ' = d r , ' = d t , and the real-valued functions q = Q(r), i = 1...4, are to be specified 
within each Szekeres class. The Szekeres metrics (with or without a cosmological 
constant) can be divided into two classes, depending on whether (Re u )' vanishes or 
not [18|, fZEfl . The class defined by (Re u )' = contains shell-crossing singularities 
which are physically "mild" — they are gravitationally weak |5U| and geodesically 



complete |JT| — and hence will not be considered here. We shall focus on the class for 
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which (Re u )' ^ 0, which admits shell-focusing curvature singularities, but no shell- 
crossings. 

The function R(t, r) obeys the evolution equation 



R 



(6) 



where M(r) is a real-valued free function. Introducing an auxiliary variable 77, defined 
by 

\fRdri = dt, (7) 
Equation (|6]) reduces to 

drj , /~6 



,w ±V^ (7±cos ^ r " 



which has solutions of the form 

"6 




r) = rio±2^—(y + l)~zF(iP,a), 

where i]q = rj(6 = 0), and F(ip,a) is an elliptic function of the third kind 
ip = sin^ 1 ty, where 



* = tanh(0/2) 



"+" in 



M> = V(7-co8he)/( 7 -l) ["-" in (I)], 



and 



where 



4 R 2 + R X - 



3M 
Ax 



X 



12M 
~^X 



sinh 9, 



7 



X 



V 

-|4/c/A|3 coshG 
-| - 4k/ sinh 6 



k > 0, 
fc < 0, 



1 



= - cosh 



-1 



3M 



-A 



A- 3 



P = \ X 2 + 



12M 



V(7-l)/(7+l)- 



(8) 

(9) 
with 

(10) 
(11) 

(12) 

(13) 
(14) 

(15) 

(16) 
(17) 



For the particular case of dust (p = 0), with gravitationally unbound matter 
distributions (k = 0), the solution reduces to the simple algebraic form: 



R{t,r) 



3M 



sinh 3 



'3A 



Mr) ~ t] 



[IB) 



Strong curvature singularities in quasispherical asymptotically de Sitter dust collapse 5 

R=0 




I 



Figure 1. Conformal diagram for asymptotically de Sitter quasispherical dust 
collapse, for the two-surface E : {x,y — const.}. The r — shell is the first one 
to become singular. The outer shell, which defines the boundary of E, first collapses 
through the past cosmological horizon (H~), and then, after a finite amount of time 
t, through the black hole event horizon (EH). The example depicted here contains a 
locally naked singularity, whose null portion is visible to local observers inside the EH. 
Since the EH coincides with the Cauchy horizon, the singularity cannot be globally 
naked. 



where t c is a real-valued arbitrary function, fixed by the initial data via 

2 



*c(0 



3A 



sinh 



-i 



Ar 3 
3M 



where the scaling R(0, r) = r was adopted. The relevant derivatives of R(t, r) are: 



R'(t, r) = R 



R(t,r) 



3M 
A 



^t' c C(.Th(/ ( -/) 



i?coth(t c — t), 



(19) 

(20) 
(21) 



where the minus sign corresponds to implosion. 

A complete solution is given by equation ([|) , where R and v are given by equations 
(!)"(© an d (jl8|) , together with the energy density 

M' - 3Mis' 



p = {M' - 3Mv')e- a - 2 P 



-2u 



(22) 



R?{R' + Rv') 

Imposing the weak energy condition f3~3fl , T a bU a u b > =3- p > 0, leads to the constraint 
1M' 
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2.1. Regularity conditions 

We assume the following regularity conditions for the spacetime: 

(i) There are no shell crossing singularities: (Re u )' > 0, which also guarantees 
that the proper area radius of the shells is a monotonically increasing function of the 
comoving spatial coordinate r, 

R' + Rv' > 0. (24) 

(ii) Regularity of the metric at r = requires that all the c$(r) be at least C 1 
therein, with finite first derivatives, such that 

lim c'Ar) < oo =>■ lim v' < oo. (25) 

r-»0+ r->0+ 

(iii) In order for the metric to be locally Euclidean at r = 0, we must have 

fc(0) = 0. (26) 

(iv) From a physical point of view, it is reasonable to require that the initial velocity 
profile, R(0, r), and the energy density, p(0,r,x,y), be everywhere bounded. With the 
scaling R(0, r) = r, it follows from equation @ that M > 0{r) near r = 0, and equation 
( fl2|) implies M' ~ 0(r 2 ); near r = we must therefore have 

M' 

~ Oir- 1 ). (27) 
3. Existence and visibility 

Each two-surface E : {t,r = const.} is a two-sphere with proper area radius 
R(t,r)e u ( r ' x ' y \ and coordinate center 2a) on E. Since the center is r-dependent, 
the shells are not in general concentric, which accounts for the absence of spherical 
(or otherwise) symmetry. Note also that, unlike in the spherically symmetric case, the 
energy density is not constant over each two-sphere. If three of the Cj's are constant, 
there is a single free function — which fixes k(r), via equation (|5|) — and the model 
reduces to the TBdS case (where the k = case corresponds to the choice c\ = c 4 = \ 
and C2 = C3 = 0). 

From equation (P2"|), one sees that the energy density diverges at R = 0, thereby 
signaling a curvature singularity As in the spherical case, r > events for which 



R(r,t) = 0, are spacelike and hence cannot be naked [19|, [2(| |34|. Of potential interest 
is the "central" r = singularity. 

In order for this singularity to be at least locally naked, there must exist at least 
one outgoing non-spacelike geodesic with past endpoint at the singularity. We note that 
the existence of a solution of the, e.g., outgoing radial null geodesies equation, does 
not per se guarantee the local visibility of the singularity. If such geodesies are emitted 
after the apparent horizon (AH) forms, they will be contained in the causal past of a 
trapped region (hence unavoidably trapped), and thus any initially diverging geodesic 
congruence will immediately start reconverging upon emission. Therefore, although 
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there may be outgoing non-spacelike geodesies with past endpoint at the singularity, if 
they are emitted after the AH forms, they are only defined at the point of emission, and 
hence we shall not consider them in the operational definition of local nakedness adopted 
in this paper. We shall show that there are outgoing non-spacelike geodesies that form 
before or at the time at which the AH does, by analyzing the geodesic behaviour along 
null and timelike directions. 



3.1. Null geodesies 

The equation for outgoing radial null geodesies (ORNGs) is 

^f] =R' + Rv'. (28) 

rfr /out 

If it admits a regular solution, t(r), in an open neighbourhood containing r = 0, the 
singularity is visible — at least locally — provided t(0) < £ah(0), where t^^r) is the AH 
curve on the {t, r} plane. Let us then assume that such a solution exists, and is given 
to leading order by 

t( r ) = t + ar a , (29) 

where a G M + , and a G N + . Expanding the central energy density near r = as 

+oo 

Pc(r)=J>^ (30) 

i=0 



we have, from equations fl22|), @, (g7|), flU) and (§), 

M(r) = M r 3 + M n r n+3 + 0(r n+4 ), (31) 

t c {r) =t + t n r n + O(r n+1 ), (32) 

R(t,r) = 3 (m*t + m n r n+1 ) (t + t n r n - t)§ + 0(r n+2+CT ), (33) 

where t n and m n are real coefficients linear in M n = (4n/n)p n , with n > 0; 
p n = (d n p c /dr n ) r=0 is the first non- vanishing derivative of the central energy density 
distribution, and 



( » = 8i,lh "(v4j- (34) 

Since the geodesic must lie on the spacetime, from equations (^) and ( |32"1) it follows 
that a > n. If o = n, we have the additional constraint, a < t n . 

The existence of a self-consistent solution of equation ( ^8|) is dependent on the 
r-differentiability of R and u at r = 0. From equation (|]), we have 

7 x c / 1 (x 2 + y 2 ) + 2(c / 2 a: + c / 3y) + ^ 

z/ (r,x,y ) = ^ . 35 

ci(x 2 + y 2 ) + 2(c 2 x + c 3 y) + c 4 

Since all the c»(r) are assumed to be at least C 1 , to leading order their MacLaurin 
expansion is 

Ci (r) = 4 + c n Sr n \ (36) 
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where c"* = (d ni Ci/dr ni ) r=Q is the first non-vanishing derivative of Cj(r), and > 1. 
Without loss of generality, we can take rii = minjrij, i = 1...4}, such that near r = we 
have 

v> = (37) 

C(x, y) = 4 + (2* + 2c°y + eg) c7 ni (x 2 + y 2 )" 1 . (38) 

If £ = (which is a set of measure zero in the initial data, since it requires c° = 0, 
% = 1...4), then v' = n\r~ x and the metric is just C° at r = 0. In general, ( ^ 0, which 
renders the metric C 1 at r = 0, as desired. In this case, from equations (p8|)-(|29D and 
([33]), we obtain, near r = 

aar^ 1 = ( — + lj + - ni R~ +ni + 0(r CT+1 -3), (39) 

9MnV i 



, , a- (40) 

Since ni > 1, the second term on the right-hand-side is to be neglected to leading order 
in r, and a self-consistent solution exists provided 

«(t + 1 J' (41) 

In 

-=l + f- (42) 

The condition a > n now reads n < 3. For n = 1,2 (i.e., for p x 7^ 0, or p x = and 
p 2 7^ 0) there is a self-consistent solution to the ORNG equation in the limit t — > t , 
r — > 0, and thus there is at least one ORNG starting from the singularity. 

Let us now examine the case n = cr = 3. Proceeding as before, we obtain, to leading 
order, 

1 



3ar 2 = 3 (^r (t3 _ a) § r 2 ; (43) 



9M n ^ 1 

T 

which is identically satisfied provided 

a 3 - p 3 a 2 + 2p 2 £ 3 a = 0, (44) 

where p = (9M0/4) 1 / 3 . This equation has two non-zero distinct roots (other than the 

4 

a = trivial root), given by a = (p 2 /2) ± y'p 4 — 8i 3 , if i 3 < | (^j 2 -) 3 , which imposes a 
constraint on p 3 , for a given p . In addition, one must also require a < t 3 , which leads 
to 

- 4 + vV + 16 - 4/i 2 < t 3 < ip 4 . (45) 

2 o 
Hence, as long as one restricts ourselves to initial data that satisfies the above condition, 
there exists an ORNG with past endpoint at the singularity. We note, however, that 
the n = 3 case is less generic than the n < 3 cases, as it requires that p\ = pi = and 
p 3 obey condition 
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Summarizing, for a given initial density profile, the ORNG equation admits a 
regular solution at r = 0, independently of the details of u(r, x, y) (provided v obeys 
the regularity condition (|25|)). 

3.2. Timelike geodesies 

Let us consider radial timelike geodesies (RTGs) described by a tangent vector field 
K a = 2§-, where r is an affine parameter along the geodesic. A sufficient set of equations 
for K a is 



K l = ±^l + e 2a (K r ) 2 , (46) 

an+^AKjg + r*±^ + <^! fl,, + fiV + fi -" = o, (47) 

V A * ; R' + Ru' K l R' + Ru' ' V ; 

where the first equation is the unit-norm condition, and the second one follows from the 
geodesic equation. By inspection, one sees that the above set admits the trivial solution 

K l = ± 1, (48) 

K r = 0, (49) 

which leads to 

t = t ±(r-T ), (50) 

r = r = const., (51) 

where r is the proper time at which RTG's depart (arrive at) the central singularity, and 
the plus or minus sign refers to outgoing or ingoing RTGs, respectively. The outgoing 
RTG departing from the singularity is given by r = 0, and t = to + r — r , and thus 
does not belong in the spacetime. The ingoing RTG is given by r = 0, t = t — r + r , 
where to = t c (0) = is the time at which the RTG arrives at the singularity. 

Since equation ( |47D is a mixed first-order linear PDE for K r (t,r), its solution need 
not be unique J35[. Indeed, one can explicitly construct other families of solutions, as 



follows. Near the singularity, we can write, to leading order 

tRTc(r) = to + br p , (52) 

fl(Wr),r) = a r 2 ~i +1 + 0{r*^ ), (53) 

#(Wr), r) = a (— + l)r 2n " + 0(^ +1 -f ), (54) 

where ao G M + , b G R, and p > n (such that the geodesies thus constructed belong in 
the spacetime). 

Let us now assume that K r (t,r) oc (t — to) a r^ ', where a,/3 G R. From equation 
(|52|), along the RTG, we have then 

K r (t RTG (r),r) = kr ap+l3 . (55) 

where G R is a constant, yet to be determined. With these ansatze, the RTG equation 
reads 

= bpr p - 1 = ^ = y/(K r )~ 2 + (R' + Ru') 2 . (56) 



Strong curvature singularities in quasispherical asymptotically de Sitter dust collapselO 



Now, from equation (|37|), we have v 1 = rii( l r ni l , where n\ > 1 is fixed by the initial 
data, and characterizes the differentiability of the metric at the origin. From equation 
(0), it follows that 

& y r 2(p-D = fc -V- 2 (^) + a^ 2 r 2 '+a^r 2 ( ,?+ " 1 - 1 )+2a^^r^ + " 1 - 1 , (57) 

where q — 1 + 2n/3, and Z = 2n/3. Now, since ni > 1, a straightforward calculation 
shows that if n < 6, the last two terms are of higher order than the remaining ones. 
We shall assume this for now, and conclude below that a self-consistent solution always 
satisfies this condition. We then look for solutions which are homogeneous in the leading 
order of r. This leads to a coupled algebraic system for the parameters p, q,a, j3: 

p - 1 = - ap - p, (58) 
ap + Q = -q + 1, (59) 

which is solved by 

2n 1 - 8 

P= 1 + ^ = TT^- (60) 
o 1 + a 

The constraint p > n now reads n < 3, which is consistent with the earlier assumption 



n < 6. For parameter values satisfying equation (]60|) , equation (|56D becomes 

C(a ,&,p,fc)r 4 "/ 3 = 0, (61) 
whose solution is given by the algebraic constraint 

4-T7 2 

C(a ,b,p,k) = b 2 p 2 -k- 2 -a 2 = 0. (62) 

9 

Now, since K r is obtained from the derivatives of R(t,r), the parameter k is not 
independent of a ; therefore, for given initial data, the set of parameters {a ,p,k} is 
uniquely determined. From equation (|62|), it then follows that b is entirely determined 
from the initial data, modulo its sign, i.e., for a a given set of initial data, both ingoing 
(6 < 0) and outgoing {b > 0) radial timelike geodesies exist. 



3.3. Local visibility 

Thus far we have shown the existence of outgoing null and timelike geodesies with 
past endpoint at the singularity. In order for these geodesies to be visible to local 
non-spacelike observers, they must not be contained in the domain of dependence of 
any trapped surfaces. Clearly, a necessary and sufficient condition for the singularity 
to be locally naked is that future-directed geodesies departing from the singularity, do 
so before or at the time at which the AH forms; otherwise they will be unavoidably 
trapped. 

The AH is the outer boundary of a trapped surface — a compact spacelike two- 
surface whose outgoing and ingoing null geodesic congruences have vanishing expansion. 
Let us then consider the two-surface £ : {t, r = const}, and a congruence of radial null 
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geodesies, with tangent vector field £ a = ^ (where A is an affine parameter along the 
geodesic), orthogonal to it. The expansion of such a congruence is given by the scalar 

e = v a c. (63) 

Orthogonality to X allows the choice 

r = e = o, (64) 

and the null-norm condition gives 

{ey = {r' + My {c? . (65) 

One can choose the affine parameter such that 

e = (R' + Rv')\ (66) 
C = e = ±1, (67) 

where the plus or minus sign corresponds to outgoing or ingoing geodesies, respectively. 
From equations (0)-© an d (|66"1)-(|6"7|), we have then 



6 = ^(R' + Ry') (R + e). (68) 



Now, from equation ( p4|) the first term is always positive-definite, and thus sign O = 
sign (R + e). Since we are interested in implosion situations, we must also have R < for 
all times. It then follows that for ingoing null geodesies (e = — 1), we have O < 0, i.e., 
their expansion is always negative (convergence). For outgoing geodesies (e = +1) one 
can have initial divergence, if R(0, r) > — 1. However, since(D -R(0, r) = — + |r < 0, 
there is a finite time £m(?") > such that R(tyi,r) < —1. Therefore, there is a time 
£ah £ (0,£m), such that 

R(t AR ,r) = -l. (69) 

It then follows from equation @ that the AH is given by the curve £ah(?~), where 
R(t\n, r) is a solution of 

—R 3 — R + 2M = 0. (70) 
3 

This equation has three distinct real roots if 3Ma/A < 1, two of which are positive and 
given by 

2 



R± = —p= sin 



^sin-^MVAH^lii; 



(71) 



where 5(x) is the Dirac delta-function, and R- > R + > 0, corresponding to the choice 

< sin -1 lo < 7r/2, < u < 1. The third root, R3 = —R- — R + , is negative and hence 
unphysical. i?_ is a generalized cosmological horizon (i?_ = a/3/A, when M = 0) and 
R + the black hole apparent horizon (R + = 2M when A = 0; the apparent and event 
horizons coincide in the static case). For 3M^/X = 1, the two horizons coincide. If 

1 Near r = this condition becomes p c (0) > A, which is always satisfied for physically reasonable 
values of the central density and cosmological constant (in geometrized units). 
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3Ma/A > 1, there is one negative real root and two complex (conjugate) roots, all of 
which are unphysical: the spacetime does not admit any horizons in this case. From 
the "+" solution in equation (|7T|), together with equations (|T8|)-(|T9|), we obtain 

tAn(r) = t c (r)-s(r), (72) 

where 



q(r) = sinh 




3A \3MVA 



1/3 



sin f— sin \3MVX) 



3/2' 



(73) 



At the origin, ^(0) = 0, and thus £ah(0) = t c (0). Since the AH and the singularity 
curve form at the same time at r = 0, the visibility of the singularity is determined by 
the relative slopes of the curves tAH( r ) and toRc( r ) on the {t,r} plane (where 'ORG' 
denotes non-spacelike outgoing radial geodesies). The singularity is (at least) locally 
naked iff 

£s(^)/(^)> 1 - ^ 

For ORNG's, given by n < 3, from equations (|29| ) and (|^), we get 

lim ( I ( = lim ^r-^n/z = + ^ (?5) 
r-*o \ dr J \ dr J r->o aa 

and hence the singularity is locally naked. This is in agreement with the analysis of f2|, 
which shows that the value of A does not alter the visibility of the singularity for n < 3. 
(For the highly non-generic case n = 3, part — but not all — of the singularity spectrum 
may be covered for sufficiently large (but finite) A; the method discussed in the present 
paper is not appropriate for such an analysis, and we thus refer to reader to reference 
@ for further details). 

For ORTG's, given by equation (]56|) , the local visibility condition reads 



nt 

\im-Pr n - p >l. (76) 

r^O bp 

Clearly, we must have p > a > n, which yields lim r ^o z ^ Lr ™ _p = +oo. Hence, the 
singularity is locally visible along such ORTG's. 



4. Curvature strength 

Let 7(t) be a RTG, with tangent vector K a = that terminates at the singularity at 
r = tq. Such a singularity is said to be Tipler strong along j(t) if the volume three- form 
V(t) vanishes in the limit r — > tq [j3~6| . If the scalar ^ = R a bK a K b obeys the strong 
limiting focusing condition, 

lim (r - r ) 2 ^ > 0, (77) 



then the singularity is gravitationally strong in the sense of Tipler ||37|| . This sufficient 
condition guarantees that any three-form defined along j(t) vanishes in the limit r — > r , 
due to unbounded curvature growth. 
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Figure 2. Locally naked singularity in quasispherical asymptotically de Sitter dust 
collapse. The dashed line is the apparent horizon (AH), which joins the event horizon 
(EH) when the boundary of the two-surface £ crosses the latter, at event P. Also 
shown is the worldline of a timelike observer O, who penetrates the EH and terminates 
at the spacelike portion of the singularity. PNC denotes the past null cone of O, which 
intersects the EH and the null portion of the singularity: clearly, there are non-spacelike 
geodesies emitted from the singularity, before the AH forms, which are contained in 
the PNC of O; the singularity is therefore locally naked. 



From equations (| 

(1- 
R 



and 



o- 1 



/ M' l.U 



V R 2 R' R 3 



P9|), we have 
A, 



R> 



(78) 
(79) 



Now, from equation (^) it follows that lim r ^ + = 0, and from equation (|23| ) we have 
lim r .^ + v' < oo. Therefore, lim r ^ + £ — 0, which implies 



lim ^ 

r-*0+ 



Using equations (^T 



M' 

lim — A. 

r^0+ R 2 R' 

3), together with equation 



we obtain 



lim * 



77 3 I f? 3 H — rat n r i] 3 



77 = t rt r" T- (r 

Hence, at r = we have 



3 
to). 



lim (r 

T— >TQ 



(80) 

(81) 
(82) 

(83) 



and the singularity is therefore Tipler strong. That is, along the timelike geodesies given 
by equations (f48"l) -([49"l), the singularity is Tipler strong, independently of the details 
of the initial density distribution. We note that, since all the sufficient criteria for 
Tipler strong singularities are sufficient for deformationally strong ones ||38[| , the central 
curvature singularity in TBdS collapse is also deformationally strong. 
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5. Concluding remarks 

We found that there are initial data, for which null and timelike families of outgoing 
radial geodesies with past endpoint at the singularity exist, which are emitted along 
the apparent horizon curve; the singularity is therefore at least locally naked. When 
examined along timelike geodesies terminating at the singularity, the latter was found to 
be Tipler strong, regardless of the initial density distribution. This constitutes a rather 
robust result, in that it holds true independently of the initial data [the initial density 
profile p(0,r,x,y), for any given velocity profile R(0, r) < 0], thereby implying stability 
against perturbations of the latter. 



It was shown in that there is no directional dependence of the local visibility 
in asymptotically flat Szekeres dust collapse. We obtained the same result for 
asymptotically de Sitter quasispherical dust collapse, as expected. We also showed 
the absence of effects of A on the visibility of singularities for n < 3, already observed 
in spherical Tolman-Bondi-de Sitter collapse. The fact that our results are qualitatively 
equivalent to those for spherical inhomogeneous collapse (with or without a cosmological 
constant), suggests that neither the lack of asymptotic flatness, nor mild (in the well- 
defined geometric sense of Szekeres metrics) departures from spherical symmetry change 
the standard picture of singularity formation and structure in gravitational collapse. 

Regarding the generality of the present results insofar as the matter content 
is concerned, two points are worth mentioning. Firstlywe considered the special 
k(r) — case (cf. equation (|])), which corresponds to gravitationaly unbound matter 
configurations. This simplifying assumption (k = 0) was found not to qualitatively 
change the corresponding results for spherical LTB collapse |39| (note, however, that 
in the asymptotically flat case, k = corresponds to marginally bound matter 
configurations). Since the structure of the central singularity in the present case is 
analogous to that of LTB collapse (the dependence of the visibility of the singularity 
on the initial data is exactly the same), we do not expect the inclusion of k{r) ^ 
to qualitatively change the end result of collapse. Secondly, a more general issue — that 
concerns not only the present analysis, but the issue of generic gravitational collapse — is 
that of the physical reasonability of "dust" . Whilst early studies of high-density nuclear 
matter suggested that such an approximation could be a legitimate one (at least for 
spherical collapse) [EOj, it is now well-known that radial and tangential pressures must 



be taken into account at the very late stages of realistic astrophysical collapse |41 
The inclusion of radial pressure in spherical dust collapse has been recently considered 
|^|, and it has been found to cover part — but not all — of the singularity spectrum, 



i.e., configurations that would otherwise develop locally naked singularities, end up 
in a black hole. The inclusion of tangential pressure has been considered for special 
cases of the Einstein cluster class [O, H3| and for more generic configurations [03 , with 



markedly different results: tangential stresses tend to uncover part of the singularity 
spectrum, i.e., configurations that would otherwise terminate in a black hole, develop a 
locally naked singularity. For the reasons mentioned above, we expect the inclusion of 
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radial and tangential stresses in quasispherical asymptotically de Sitter dust collapse to 
produce similar effects to those of the spherical case. 

Unlike the matter content, however, the effects of non-sphericity on the endstate of 
collapse are much less clear, and further analyses, of stronger deviations from spherical 
symmetry, need to be undertaken to confidently establish the role of the latter in the 
final state of generic gravitational collapse. 
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